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A  Comparison  of  Test  Methods  for  Determination 
of  Flexural  Strength  in  Urea  Model  Ice 

SHARON  L.  BORLAND 


INTRODUCTION 

Studies  of  forces  and  displacements  resulting  from 
the  interaction  of  ice  and  structures  commonly  involve 
physical  models.  The  accuracy  of  such  models  to  simu¬ 
late  ice-structure  interaction  relies  heavily  on  their  abil¬ 
ity  to  correctly  model  both  geometry  and  mechanical 
properties.  In  many  practical  problems  involving 
ice-structure  interaction  (e.g.  artificial  islands  for  petro¬ 
leum  exploration  in  the  Arctic,  conical-shaped  light 
piers,  ice-breaking  ships),  bending  or  flexure  of  the  ice 
is  an  important  failure  mode.  This  makes  the  flexural 
strength  a(  one  of  the  most  significant  properties  of  the 
modeling  material.  For  example,  in  the  CRREL  test 
basin  when  urea-doped  ice  is  used,  the  ice’s  flexural 
strength  is  monitored  throughout  the  tempering  period 
of  an  ice  sheet  (approximately  20  hours).  Testing  pro¬ 
ceeds  only  when  the  desired  value  of  a(  is  reached. 
Because  of  the  strong  time  dependency  of  orfor  urea  ice, 
the  method  for  measuring  <jf  must  be  quick,  accurate  and 
efficient. 

An  in-situ  cantilever  beam  test  (CBI)  method  is 
currently  used  to  determine  the  flexural  strength  of  ice  in 
the  CRREL  test  basin.  This  test  measures  the  load  P 

max 

necessary  to  break  the  beam.  Flexural  strength  is  then 
calculated  usinsi  P  and  elastic  beam  theory.  The  CBI 
test,  used  extensively  to  measure  of  for  ice  in  the  field  as 
well  as  in  the  laboratory,  is  relatively  simple  to  set  upand 
perform,  and  it  provides  consistent  results  when  done 
properly.  This  method,  however,  yields  0(.  values  that 
may  be  considered,  at  best,  to  be  an  index  value  of  true 
flexural  strength  of  the  ice.  The  actual  flexural  strength 
can  only  be  calculated  by  recognizing  an  ice  beam  as  a 
nonhomogeneous,  anisotropic,  elasto-viscoplastic  ma¬ 
terial.  Unfortunately  this  is  impossible  with  current 
analytical  techniques.  Good  approximations  of  afcan  be 
made  using  elastic  beam  theory.  Also,  certain  con¬ 
straints  in  testing  technique  (as  outlined  in  the  recom¬ 
mendations  of  the  I AHR  Working  Group  on  Ice  Testing 
Methods  (Schwarz  et  al.  1981  )|  must  be  followed. 


Another  method  used  to  determine  ofis  the  simple 
beam  test  with  three-point  loading  (SB).  This  test  can  be 
conducted  in  situ  (SBI)  or  out  of  the  water  (SBO). 
Problems  with  beam  root  deformations  due  to  plate  ef¬ 
fects  and  stress  concentrations  in  CBI  tests  can  be 
minimized  with  SB  tests  (Svec  and  Frederking  1981). 
Furthermore,  in  at  least  one  case  [tests  on  urea  ice  in  the 
NRCC  ice  model  tank  (Timco  1985)),  results  from  SB 
tests  have  been  found  to  be  the  same  as  those  from  CBI 
tests.  Since  the  SB  test  may  be  more  reliable  for 
determining  a(..  I  performed  a  series  of  experiments  to 
compare  0r  obtained  with  CBI,  SBI  and  SBO  test  meth¬ 
ods.  These  tests  were  conducted  on  urea  ice  in  the 
CRREL  test  basin.  Specifically  the  following  objectives 
were  defined:  1 )  determine  the  flexural  strength  of  urea 
ice  using  the  two  beam  methods,  2)  compare  the  three- 
point  simple  beam  method  for  both  in-situ  and  out-of- 
waler  conditions,  and  3)  compare  a(.and  its  variation 
with  respect  to  beam  geometry. 

EXPERIMENTAL  PROCEDURE 

Ice  sheet  growth  and  description 

The  urea  ice  sheets  used  to  produce  the  beam  samples 
in  this  study  were  grown  from  an  aqueous  solution  of 
0.95%  urea  by  weight.  To  assure  uniform  distribution  of 
thickness  and  mechanical  properties  throughout  each  ice 
sheet,  the  following  techniques  were  used  to  produce 
each  ice  sheet  tested.  The  solution  was  initially  mixed 
using  airbubblers  and  circulating  pumps  and  then  cooled 
by  heat  exchange  at  the  air/water  interface  until  the 
solution  reached  a  uniform  temperature  of  -0. 1  °C.  The 
bubblers  and  pumps  were  then  shut  off,  and  the  solution 
was  "wet  seeded"  by  spraying  a  fine  mist  of  water  into 
the  cold  (- 1 2°C)  air  above  the  water  surface.  The  result¬ 
ing  ice  crystals  settled  on  the  water  surface  and  nucleated 
the  ice  sheet.  The  ice  sheet  was  grown  at  an  ambient  air 
temperature  of  -I8°C  until  a  target  thickness  was 
achieved.  At  this  time  the  air  temperature  was  raised  to 


Figure  I .  Test  set-up  and  corresponding  beam  diagrams. 


-4°C  to  retard  ice  growth  and  to  help  stabilize  the  ice 
properties.  A  more  detailed  discussion  of  ice  growth 
techniques  for  the  CRREL  test  basin  is  provided  by 
Hirayama(1983). 

The  resulting  ice  sheet  had  a  two-layered,  columnar 
structure  similar  to  most  ice  sheets  grown  using  this 
technique  (Hirayama  1983,  Gow  1984,  Tirnco  1984). 
Typically  the  top  layerconsisted  of  fine-grained  colum- 
narcrystals  and  was  thinner  and  stronger  than  the  bottom 
layer.  The  bottom  layer  was  made  up  of  coarse-grained 
columnar  crystals.  Ice  crystals  ranged  in  sLe  from  1  to  2 
mm  in  the  top  layer  and  were  5  mm  or  more  in  the  bottom 
layer  for  a  typical  7-cm-thick  ice  sheet. 

Beam  test  methods 

A  total  of  23  series  of  beam  tests  were  systematically 
conducted  on  iOdifferent  ice  sheets  ranging  in  thickness 
from  about  35  mm  to  1 10  mm.  Each  series  consisted  of 
four  consecutive  beam  tests  and  took  about  an  hour  to 
complete.  The  first  test  in  a  series  was  the  CBI.  followed 
by  the  SBI  and  the  SBO.  An  additional  CBI  test  was  done 
as  the  fourth  test  in  the  series  to  document  any  change  in 
the  ice  strength  over  the  time  span  of  testing.  Beams  were 
failed  in  top-tension  (i.e.  the  top  surface  of  the  ice  sheet 
was  in  tension)  and  bottom-tension  modes.  Initially  the 
load  actuator  speed  was  varied  from  6  to  1 2  mm/s.  but  no 
effect  on  the  strength  values  was  observed  so  the  rate  was 
kept  at  1 2  mm/s.  This  speed  resulted  in  a  tinie-to-failure 
of  0.5-1  s. 

The  CBI  tests  were  conducted  using  guidelines  rec¬ 
ommended  by  the  IAHR  Working  Group  on  Testing 
Methods  in  Ice  (Schwarz  et  al.  1981)  and  Tirnco  (1981). 
The  beam  geometry  (Fig.  1 )  was  described  in  terms  of 
the  ratio  of  length  L  to  thickness  h  (5:1  to  7:1  for  my 
beams,  as  recommended  by  Tirnco  for  model  ice)  and  the 
ratio  of  width  B  to  thickness  (1:1  to  2: 1 ).  This  geometry 
has  been  shown  to  minimize  buoyancy  effects  (Tatin- 
claux  and  Hirayama  1982).  Cantilever  beams  were  hand- 
sawn  from  the  ice  sheet  (Fig.  la).  Each  beam  was  broken 
in  downward  or  upward  bending  by  using  a  motor- 
driven  actuator  to  apply  load  P  to  the  tip  of  the  beam.  A 
load  cell  mounted  on  the  actuator  was  used  to  monitor 
the  maximum  applied  force  /*  at  failure.  After  each 
beam  was  broken.  L  and  B  were  measured  with  a  ruler  to 
an  accuracy  of  1  mm,  and  h  was  measured  with  a  vernier 
caliper  to  an  accuracy  of  0. 1  mm.  The  flexural  strength 
0fwas  then  calculated  using  simple  elastic  beam  theory 
from 

Or  =  6  _  .  ( 1 ) 

Blr 

Typically  five  beams  were  tested  in  rapid  succession  and 
the  results  averaged. 

The  SB  tests  were  conducted  using  the  same  IAHR 


guidelines  as  for  the  CBI  tests  with  but  a  few  exceptions. 
Beams  were  cut  out  of  the  ice  sheet  and  loaded  at  three 
equidistant  points  in  the  SB  test  apparatus  (Fig.  lb.c). 
After  the  beams  were  cut  out  of  the  ice  sheet  and  tested, 
I  discovered  that  about  a  third  of  the  beams  had  L/h  ratios 
less  than  5:1 . 

l  he  SB  test  apparatus,  designed  for  laboratory  tests, 
has  two  roller  supports  and  a  slot  at  the  top  to  ensure 
proper  placement  of  the  loading  point  at  the  midspan  of 
the  beam.  The  length  between  the  roller  supports  is  ad¬ 
justable  to  accommodate  various  beam  lengths.  The 
beam  lengths  were  generally  250  .  330  and  450  mm, 
depending  on  the  ice  thickness.  The  same  motor-driven 
actuator  used  for  the  CBI  tests  was  used  for  the  SB  tests. 

For  the  SBI  tests  the  apparatus  was  placed  in  the  water 
and  supported  in  a  slot  cut  in  the  surrounding  ice  sheet 
(Fig.  1  b).  Wood  blocks  were  wedged  between  the  frame 
and  the  ice  sheet  to  provide  adequate  clearance  to  just 
slide  the  floating  beam  over  the  roller  supports.  After 
aligning  the  beam  in  proper  position  in  the  frame,  a 
downward-acting  load  was  applied  to  the  top  of  the  beam 
at  the  midspan.  To  cause  bending  in  the  beams  so  that  the 
fibers  at  the  top  ice  surface  were  stressed  in  tension  (top- 
tension  tests),  the  beams  were  flipped  over  before  plac¬ 
ing  them  in  the  frame.  Similarly,  to  stress  the  bottom 
fibers  in  tension  (bottom-tension  tests),  the  beams  were 
placed  right  side  up  in  the  frame.  The  flexural  strength 
was  calculated  using  elastic  beam  theory: 

T,  P  I 

°t  =  (2) 

2  Blr 

where  P  =  maximum  force  at  failure 

max 

L  =  length  between  the  roller  supports 
B  =  ice  beam  width 
h  =  thickness  of  the  beam. 

Since  the  beam  was  tested  in  the  water,  the  effect  of  the 
weight  of  the  beam  was  cancelled  out  by  buoyancy  ef¬ 
fects.  Several  beams  were  tested  in  each  series  and  the 
results  averaged  in  each  test. 

The  SBO  tests  (Fig.  I c)  took  place  much  like  the  SBI 
tests,  the  only  difference  being  that  tests  were  accom¬ 
plished  on  a  carriage  just  above  the  water  surface.  The 
beams  were  removed  from  the  water  and  then  carefully 
placed  and  aligned  in  the  frame  as  quickly  as  possible. 
The  flexural  strength  was  then  calculated  from 

0r  =  j-  Pmax  L  +  "  #3) 

2  Blr  4/j 

where  p  is  the  ice  density  (920  kg/m1)  and  g  is  accelera¬ 
tion  due  to  gravity.  Again,  as  for  the  previous  beam  test 
methods,  several  beams  were  tested  and  their  results 
averaged. 
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Table  1.  Flexural  strength  results  from  each  test  series.  Each  value  represents  the  mean  result  from  3-5  tests, 
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Figure  2.  Flexural  strength  values  obtained  from  initial  and  final  in-situ  cantilever 
beam  tests  for  each  series.  Each  data  point  represents  the  mean  value  of  several  beam 
samples,  and  the  corresponding  standard  deviations  are  indicated  by  the  bars  at  each 
point. 


RESULTS 

Comparison  of  cantilever  beam  tests 

The  values  of  a,,  from  the  initial  and  final  cantilever 
beam  tests  in  a  test  series  (Table  I ).  for  both  top-  and 
bottom-tension  modes,  are  shown  in  Figure  2.  These 
data  indicate  good  one-to-one  correspondence  between 
the  CB1  tests.  Using  a  linear  regression  analysis,  the 
equation  of  the  best-fit  line  was  found  to  be 

Or,  =  0.83  (o,,)  +  8.4  (4) 

with  9 1  %  correlation  (;  '  =  0.82).  Although  the  slope  of 
this  line  suggests  that  the  ice  sheet  strength  was  reduced 
slightly  over  the  time  it  took  to  run  a  test  series,  this 
change  may  be  considered  negligible  since  the  variation 
coefficient  of  each  data  point  is  10-15%.  This  indicates 
that  a  slope  of  I  and  an  intercept  of  0  (i.e..  no  change  in 
strength)  can  reasonably  be  assumed.  Thus,  the  ice 
sheet's  mechanical  properties  were  essentially  stable  for 
the  duration  of  each  test  series.  Both  top-  and  bottom- 
tension  test  results  show  this  behavior. 

Comparison  between  simple  beam  tests 

A  plot  of  the  results  from  the  SBI  and  SBO  tests  (Fig. 
3)  indicates  excellent  one-to-one  correspondence  be¬ 
tween  the  two  test  methods.  Regression  analysis  re¬ 
vealed  a  correlation  coefficientof  95%  (/"  =  0.90)  for 


straight-line  fit.  It  appears,  then,  that  there  is  no  signifi¬ 
cant  difference  in  a,,  results  from  using  either  test  method. 
However,  the  o,.  values  ranged  from  20  to  100  kPa,  and 
this  conclusion  may  only  be  valid  within  this  range.  This 
conclusion  holds  true  when  the  top  surface  of  the  ice 
sheet  was  in  tension.  Analysis  of  the  few  bottom  tension 
tests,  however,  shows  a  slightly  higher  slope  (1.21)  for 
the  best-fit  line.  These  results  indicate  that  the  brine 
drainage  resulting  from  the  removal  of  the  ice  beam  does 
not  significantly  alter  the  structural  integrity  of  the  ice 
and  therefore  its  flexural  strength. 

Comparison  between  simple  beam 
and  cantilever  beam  tests 

A  comparison  of  c,.  values  from  each  simple  beam 
test  with  the  average  a,  value  front  the  initial  and  final 
cantilever  beam  tests  is  shown  in  Figure  4.  The  plotted 
values  of  a,  (for  both  top-  and  bottom-tension  modes) 
indicate  a  good  one-to-one  correlation.  Clearly  there  is 
more  scatter  in  Figure  4  than  in  Figures  2  or  3.  Most  of 
the  scatter  seen  here  is  the  result  of  combining  the  data 
from  Figures  2  and  3:  the  small  amount  of  scatter  seen  in 
these  figures  is  amplified  in  Figure  4.  Some  of  the  scatter 
may  also  lie  due  to  preslressing  of  the  beams,  which 
unavoidably  occurred  when  handling  the  fragile  ice 
beams.  This  was  especially  true  for  the  out-of-water 
beam  tests.  The  correlations  here  are  not  as  strong  as  for 
the  comparison  of  the  simple  beam  methods  ( Fig.  3):  the 
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Figure  3.  Effect  of  ice  beam  placement  on  flexural  strength  results  from  the  three- 
point  simple  beam  test  method. 


correlation  coefficient  was  found  to  be  slightly  less, 
90%.  Ice  test  results,  however,  typically  show  10-15% 
variations  even  under  very  controlled  conditions.  Thus, 
these  data  indicate  that  there  is  little  difference  in  cr 
when  obtained  with  SB  or  CBI  test  techniques  either  in 
situ  or  out  of  the  water. 

Previous  results  in  freshwater  ice  (e.g.  Gow  et  al. 
1978)  show  that  significant  stress  concentrations  at  the 


root  of  a  cantilever  beam  result  in  a  much  lower  flexural 
strength  value  than  would  be  obtained  from  simple  beam 
tests.  Conversely,  previously  published  results  for  sea 
ice  (Frankenstein  1968,  Vaudrey  1978)  show  no  differ¬ 
ence  in  flexural  strength  values  from  either  test  method. 
Schwarz  and  Weeks  (1977)  pointed  out  that,  in  sea  ice, 
the  stress  concentrations  in  the  cantilever  beam  may  be 
relieved  through  the  plastic  flow.  Because  urea  ice  (used 
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Figure  4.  Effect  of  test  method  on  flexural  strength. 

6 


Figure  5.  Effect  of  beam  lengtli-to-thickness  ratios  of  the 
simple  beam  tests  on  the  relati ve  difference  between  beam  test 
methods  for  all  the  data. 


in  this  study)  exhibits  similar  plastic  flow  behavior,  the 
effect  of  the  stress  concentrations  in  the  cantilever  beam 
on  the  flexural  strength  are  probably  minimal.  Thus  the 
flexural  strength  values  as  measured  using  the  CBI  test 
or  the  SB  test  should  be  the  same  for  urea  ice  as  con¬ 
firmed  by  my  results  (Fig.  4).  These  results  also  confirm 
similar  findings  published  by  Timco  ( 1 985)  for  the  urea 
model  ice  at  the  NRCC  ice  tank. 

Effect  of  Llli 

As  mentioned  previously,  about  one-fourth  of  the  SB 
tests  were  inadvertently  conducted  using  L/h  ratios  of 
less  than  4.  These  low  L/h  ratios  were  caused  by  acciden¬ 
tal  ly  failing  to  adjust  the  position  of  the  rol  ler  supports  on 
the  SB  apparatus  to  accommodate  thick  (//  >  80  mm) 
beam  specimens  for  six  test  series.  The  (L/h)SB  ratio 
ranged  from  approximately  2  to  7  for  the  simple  beam 
tests,  while  ( L/h)CB  ranged  from  4.5  to  7  for  CBI  tests. 

The  effect  of  (L///)SB on  the  aSB/aCB  strength  ratio  is 
illustrated  in  Figure  5.  This  plot  includes  ail  the  test 
series  for  both  in-situ  and  out-of-water  simple  beam 
tests;  each  point  represents  the  mean  value  of  several 
beam  tests  conducted  in  rapid  succession.  (The  standard 
deviations  were  omitted  from  the  plot  for  clarity.)  Since 
the  abscissa  contains  L/h  values  from  the  simple  beam 
test  results.  Figure  5  shows  the  effect  of  varying  beam 
geometry  for  the  simple  beams  only.  Linear  regression 
analysis  of  these  data  indicates  a  definite  decreasing 
trend  of  the  strength  ratio  ac„/a„D  with  increasing 
( L/Ii)sb .  even  for  L/li  in  the  5-7  range. 

The  equation  for  the  best-fit  curve  in  Figure  5  was 
determined  to  be 

°SB/°CB  =  I  -38  -  0.096  (Z.//»)sb  .  (5) 

Based  on  this  regression  equation,  then,  the  value  of 


(L///)sb  that  best  corresponds  to  a  strength  ratio  CTSB/a(  B 
of  unity  is(L///)SB  =  4.  Allowing  10%  deviation  in  either 
direction  for  the  ratio  which  is  well  within 

^  D  v-  D 1 

acceptable  limits  for  ice  testing,  one  can  expect  to  be 
reasonably  confident  that  the  ratio  oSB/a(;B  is  unity  for 
(Z.//;)SB  values  of  3-5. 

It  is  possible  that  the  effect  of  (Z.//;)SB on  the  strength 
ratio  in  Figure  5  could  be  due  to  variations  in  aC(J  with  ( L / 
//)CB.  To  confirm  that  gsb  is  the  effect  and  not  gcb.  I 
plotted  a  subset  of  the  data  where  a  narrow  range  of  (LI 
h)CR  was  present.  In  Figure  6,  data  are  plotted  using  only 
those  points  for  which  (L/Ii)cb  is  restricted  to  the  narrow 
range  of  4.5-5.5.  out  of  the  full  range  of  4.5-7.  Linear 
regression  analysis  reveals  the  equation  of  the  best-fit 
curve  to  be 

cSb/°cb  =  1-6-0.145  (L/h) SB.  (6) 

Although  the  best-fit  curve  in  Figure  6  has  a  slightly 
steeper  negative  slope  than  the  curve  in  Figure  5,  the 
general  trend  of  this  data  subset  is  much  the  same  as  that 
for  the  entire  data  set  seen  in  Figure  5.  Based  on  the  above 
regression  equation  for  the  data  subset,  the  strength  ratio 
is  unity  for  (L//i)s  =4.1,  which  is  nearly  the  same  result 
obtained  for  the  entire  data  set.  thus  confirming  that  there 
is  no  appreciable  effect  due  to  the  geometric  variations  in 
the  CBI  tests. 

The  effect  of  (L/h)co on  the  strength  ratio  between  in- 

o  D 

situ  and  oul-of- water  simple  beam  tests  was  also  checked. 
The  ratio  aSBO/°SB|was  plotted  against  (L//;)SB(Fig.  7). 
This  shows  an  average  ratio  close  to  unity,  indicating  that 
(L///)SB  affects  these  two  test  methods  identically.  Note 
that  (L///)sb  was  identical  for  the  SBO  and  SBI  tests 
within  each  series. 

These  tests  have  shown  that  o(.  may  be  more  sensitive 
to  L/h  than  previously  thought.  The  effect  of  L/h  both 
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Figure  6.  Effect  of  beam  lengtli-to-tliickness  ratios  of  the  simple  beam  tests  on 
the  relative  difference  between  beam  test  methods  for  the  data  corresponding  to 
cantilever  beam  tests  in  which  (L/hl c/j  =  4.5  to  5.5. 


(L/h)SB 


Figure  7.  Effect  of  beam  lengtli-to-tliickness  ratios  on  the 
relative  difference  in  flexural  strength  values  obtained  from 
the  in-situ  and  out-of-water  simple  beam  tests. 


inside  and  outside  the  range  5-7  needs  further  inves¬ 
tigation.  Based  on  my  analysis  of  the  effects  of  beam 
geometry  on  strength  results.  I  recommend  that 
(L///)sr=  4  be  used  in  SB  tests  to  get  the  same  results  as 
from  CB1  tests  for  urea  model  ice. 


SUMMARY  AND  CONCLUSIONS 

A  total  of  ten  urea  ice  sheets  ranging  in  thicknesses 
from  35  to  110  mm  were  tested  for  flexural  strength 
using  three  beam  test  methods.  Ice  strength  was  deter¬ 


mined  by  loading  the  ice  beams  to  failure  in  flexure  using 
various  loading  configurations.  These  configurations 
included  the  current  standard  methods  using  in-situ 
cantilever  beams  and  both  in-situ  andout-of- water  three- 
point  simple  beams.  For  flexural  strength  values  of 
between  20  and  1 00  kPa.  little  practical  difference  in  o( 
was  found  between  any  of  the  beam  test  methods.  Based 
on  laboratory  observations  and  analysis,  the  following 
conclusions  were  drawn  from  this  study: 

•  The  ice  sheet  remained  stable  with  respect  to  C( 
over  the  hour  or  so  it  took  to  conduct  a  test  series, 
when  the  ambient  temperature  was  — l°C. 
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•  No  difference  in  a  .in  the  range  of  20- 1  (K)  kPa  was 
observed  for  the  three-point  loaded  simple  beam 
tests,  whether  the  beam  was  tested  in  situ  or  out  of 
the  water  using  a  three-point  beam  loading  appa¬ 
ratus. 

•  Little  practical  difference  was  found  between  three- 
point  simple  beam  test  results  and  in-situ  cantile¬ 
ver  beam  test  results.  This  is  in  agreement  with 
previously  published  results  (Timco  1985). 

•  Values  for  a,  from  the  three-point  simple  beam 
tests  were  observed  to  increase  w  ith  decreasing 
L/h  for  ratios  between  2  and  7. 

•  For  urea  model  ice.  a  beam  geometry  of  L/h  -  4 
should  be  used  forthree-point  loaded  simple  beam 
tests  to  obtain  the  same  results  as  for  in-situ 
cantilever  beam  tests. 
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